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of the bimodal distribution reveals the existence of a fixed point that is connected to itself by a homoclinic 
trajectory. The mathematical conditions associated with bimodality have been established. The critical factors 
necessary for both a breaking of symmetry and a transition from bimodal to unimodal behaviour, in the manner 
of a bifurcation, have been analysed. 

PACS numbers: 47.20.Ky, 45.70.Vn, 05.45.-a, 02.30.Oz 

Keywords: Bifurcation and Symmetry Breaking, Traffic Flow, Nonlinear Dynamics 



I. INTRODUCTION 

Bimodal distributions are characterised by the existence of 
two distinct modes (peaks) in a single frequency distribu- 
tion (ifl. Usually such a distribution is seen as a mixture of 
two normal distributions. Many natural phenomena, to think 
of just a few cases |2|-|4|], exhibit bimodality. This work is an 
attempt to provide a global mathematical model for bimodal- 
ity in a man-made context — vehicular traffic flows. 

Any standard approach to understanding bimodal distribu- 
tions refers in essence to statistical methods, and it is only 
in a rare few instances that efforts have been made to pro- 
vide a model description for bimodality [5]. In keeping with 
the latter principle, the primary emphasis of this work is on 
developing a mathematical model to furnish a single global 
description of bimodal phenomena. The model equation thus 
obtained, is then applied to comprehend the dynamics under- 
lying a bimodal distribution, following standard methods of 
nonlinear dynamics laLZl]. In so doing, the dependence of bi- 
modal properties on various parameters has been brought to 
the fore. The dynamics of the model bimodal distribution re- 
veals that the controlling parameters can be adjusted to cause 
bimodal-to-unimodal transition (akin to a bifurcation) and a 
breaking of symmetry in the system (evidence for both of 
which has been adduced in some previous studies on bimodal- 
ity Hi). 



II. BIMODALITY IN VEHICULAR TRAFFIC FLOWS 

Nowadays the number of vehicles on the roads has in- 
creased prodigiously, along with an accompanying necessity 
to expand road networks. Any implementation of an advanced 
transportation system, therefore, needs a clear view of the 
present intricacies in the fundamental patterns of traffic flows. 
Studies on vehicular traffic flows are usually aimed at un- 
derstanding interactions among streams of vehicles (modelled 
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occasionally as particles), and a substantial body of relevant 
work, ranging across multiple perspectives, has already come 
into existence Ill0l428ll . The practical motive behind these 
studies is to eliminate the problem of traffic congestion, and 
to devise effective methods of controlling traffic flow. 

The mathematical model that has been forwarded in this 
work, takes a macroscopic view of the vehicular traffic of a 
city at a given hour in the day — in effect, the bulk of traffic 
as a function of time. The data needed for the purpose of mod- 
elling have been taken from the repository of the Alabama De- 
partment of Transportation (ALDOT) and the particular city 
whose traffic data have been used here, is Jackson in Alabama 
State, USA. The principal factor guiding the choice of this city 
is that with Jackson being comparatively small in size, the data 
pertaining to it were systematically recorded, and free of the 
usual complications associated with the traffic of very large 
cities. This has the advantage of simplicity, where placing the 
foundations of a mathematical model is concerned. 

The Jackson data indicate that the bulk traffic is dominated 
by flows in the eastward and westward directions. The plot of 
the data for the westward traffic is shown in Fig. Q] with the 
broken line. Similarly, the plot of the eastward traffic is shown 
in Fig. |2 In both these plots, the time, i, at which the traffic 
volume, N (i), was measured, has been scaled so that t = 
coincides with midday, 12:00 hours. This enables one to take 
advantage of any possible symmetry in the distribution of the 
data, especially if that symmetry were to exist about t = 0. 
From a look at both Figs.Q]&[2] however, only an asymmetric 
bimodal distribution is to be seen, with two peaks of unequal 
height in either plot. What Fig. Q] suggests is that the morning 
traffic volume is much higher in the western direction than 
in the opposite direction, while the behaviour of the bimodal 
curve in Fig.|2]is quite the opposite, with the eastward traffic 
being disproportionately high in the evening. One may then 
surmise that the working quarters are in the western side of 
the city, and the residential quarters are in the east. 

It is also important here to give some careful thought to the 
way extraneous traffic (due to vehicles coming from outside 
the city, passing through it, and leaving it ultimately) may add 
to the count of the native traffic flow. To address this mat- 
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ter, it is expedient to imagine Jackson city as a bounded vol- 
ume, through which there occurs a bulk flow of traffic from 
the rest of the world. This flow is bi-directional, and along 
either direction, it is conceivable that the external traffic vol- 
ume maintains a uniform average value over time. Besides, 
whatever goes into this idealised bounded system at one end, 
also comes out of the other end in finite time. Consequently, 
the alien traffic does not cause any appreciable bias in the bi- 
modal distribution. The variation in the data around this con- 
stant background of global traffic, is then due only to the local 
city traffic, which oscillates back and forth between the east- 
ern and western boundaries, once in a cycle of twenty four 
hours. 

It will be quite in order also to mention that social and eco- 
nomic conditions have a major role to play in the way the data 
are distributed from one day of the week to the other. One 
can easily imagine that the weekend traffic (on Saturdays and 
Sundays) will follow a different pattern, in comparison to the 
traffic on a working day. And even on working days close to 
the weekend, like Mondays and Fridays, the traffic flow may 
be different from what it would be in the middle of the week. 
So a judicious choice of the data has to be made here. For 
the purpose of crafting a mathematical model from first prin- 
ciples, the data employed in this work are the average values 
of all the Wednesday traffic data of Jackson city in the year, 
2007. Wednesday is in the middle of the week, and it is in this 
period of the week that all city traffic is expected to rid itself 
of the weekend effect and reach a local mid-week equilibrium. 
So it stands to reason that the corresponding data should con- 
vey more unambiguous trends in the traffic flow. 



III. A MATHEMATICAL MODEL FOR BIMODALITY 

The continuous curves in Figs. Q] & [2] follow the model 
function posited here to describe bimodal distributions. In 
such distributions one encounters three turning points. So it 
seems natural to suggest that a fourth-degree (quartic) polyno- 
mial [8, 9], whose derivative is of the third degree (cubic), will 
correspond to the three turning points. This particular feature 
can be captured by an "inverted Mexican hat"-type function of 
the form, N (t) = a + t 2 — bt 4 , in which TV is volume of the 
traffic flow, and t is the hour at which the traffic flow has been 
measured. The discrepancy between this function and the pat- 
tern followed by the data, however, lies in the asymptotic con- 
ditions of the latter. Both the inner and the outer time limits 
on the horizontal axis in Figs.Q~|&[2] are close to midnight. At 
that hour, on a working day in the middle of the week, the vol- 
ume of traffic trickles down to small values, a condition that 
should, therefore, be described mathematically by a function 
with an asymptotic decay, N — > 0, on the extreme ends of 
the time axis. And this is exactly how it should be, if one re- 
calls that a bimodal distribution is usually treated as a mixture 
of two normal distributions, which have a similarly decaying 
tail. This, however, is not indicated by the "inverted Mexican 
hat" function, which intersects the line N = in finite time. 

To capture the required asymptotic behaviour of the bi- 
modal distribution, therefore, it is expedient to employ a 
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FIG. 1: Bimodal distribution of traffic flow due west. The vertical 
axis measures the volume of traffic, N, for a given time, t (scaled 
in hours). The horizontal axis has been scaled so that t — is set 
at midday, 12:00 hours. The dotted broken curve joins the real data 
points, and the continuous curve traces the model given by Eq. {2}. 
The parameter values used for the fit are A = 44.0, /j, = 8.53, 
A = 0.19 and /3 = -0.09. 
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FIG. 2: Bimodal distribution of traffic flow due east. The plotting 
follows the same principles outlined in Fig.Q] but the parameter val- 
ues necessary for fitting the data are A = 44.1, \x — 10.5, A = 0.22 
and j3 = 0.24. The sign and the absolute magnitude of ft make the 
most significant difference between the two bimodal distributions. 



Maxwell-Boltzmann type of function in the form, N (£) = 
a (jjl + i 2 ) cxp (— A 2 i 2 ). From the foregoing function, for 
low values of t, on approximating cxp (— A 2 t 2 ) ~ 1 — X 2 t 2 , 
one does, in fact, recover the features of the "inverted Mexi- 
can hat" function. Another aspect of the Maxwell-Boltzmann 
distribution is that its integral (the area under the curve) is a 
constant quantity, a feature in consonance with the fact that in 
an average sense the total volume of the internal city traffic is 
arguably constant over one full day. 

With the required asymptotic condition satisfied thus, it is 
now instructive to examine the role of the parameter, fi. From 
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the Maxwell-Boltzmann function, setting its time derivative, 
N = 0, one obtains the condition, t [l - A 2 (fi + t 2 )] = 0, 
which indicates three turning points in time, at t = and 

1 /2 

£ = ± (A~ 2 — /i) Now, the two non-zero roots of t can 
only be real when /i < A~ 2 . When = A~ 2 , all the three 
roots coalesce at t — 0. Noting that bimodality implies the 
existence of three turning points, the importance of the pa- 
rameter /i in maintaining bimodality is now easily evident. A 
detailed analysis of this particular issue will be presented in 
Section HVl 

The distribution traced out by the Maxwell-Boltzmann 
function is, however, symmetric about t = 0, which is cer- 
tainly not to be expected in a real bimodal distribution. And 
more pertinently for the case of traffic flows, this lack of sym- 
metry is exchanged over a time scale of half a day, as it has 
been shown in Figs.Q]&[2] While the traffic due west is higher 
in the morning hour than at any other time of the day, the traf- 
fic due east displays an exact reversal of this trend. To account 
for these two features, two new parameters, 7 and er, are now 
introduced in the symmetric Maxwell-Boltzmann distribution, 
which, consequently, goes as 



A (t) = a (fj. + t 2 ) exp (-\ 2 t 2 ) (1 + 7 ) c 



(1) 



Any non-zero value of 7 in the foregoing expression will 
raise one peak in the bimodal distribution, and lower the other, 
i.e. break the symmetry about t = in Figs.Q]&[2] Symmetry 
is, of course, restored when 7 = 0, and in this special case 
it may be noted in passing that Eq. ((TJ assumes the general 
mathematical form of the energy eigenfunction of the second 
state of the linear harmonic oscillator [29!]. When 7 / 0, 
the part played by a in the distribution given by Eq. ((TJ is 
to exchange the respective positions of the two uneven peaks, 
depending on the direction of the traffic flow, a condition that 
can be achieved simply by assigning the values, a = ±1. 

With all the parameter values specified suitably, Eq. ([TJ can 
be used to fit the data, but an expedient transformation of it 
helps in reducing the number of parameters in the model. Ac- 
cordingly, Eq. (HJ is recast as 



N(t) =A(fi + t 2 ) exp - (A* - /3) 



(2) 



in which /? = a (2A) _1 In (1 + 7) and A = a exp (j3 2 ). Now 
P, having taken over the roles of 7 and a, controls both the 
breaking of symmetry and the exchange of the peaks in the 
bimodal curve. The former condition is obtained from the 
absolute magnitude of f3, while the latter from its sign. These 
values, as well as values of A, /i and A are used to fit the data, 
as shown in Figs. |TJ&[2] 



IV. DYNAMICS OF BIMODALITY 

The model for the data fit, as given by Eq. (fJJ, is in the form, 
N = N (t). It has been realised by now that two parameters 
in this model, [i and f3, play a crucial part in causing an asym- 
metric bimodality. And it is also clear that for certain values 



of these parameters (for instance [i = A -2 and (3 = 0), the bi- 
modal distribution undergoes qualitative changes. It is worth 
investigating these aspects of the model in greater detail. 

To that end, it is necessary to find N, and then to plot 
a phase portrait, N—N. This will help in attaining an un- 
derstanding of the stability, creation and annihilation of the 
turning points, and also the properties of any possible fixed 
point(s) in the bimodal system [6, 7]. So going back to Eq. (ffj), 
its first derivative in time (indicated by an overdot) is obtained 
in a non-autonomous form as 



A 



t 



f 2 



A (At - 0) 



N . 



(3) 



Turning points are given by the condition N = 0. In ad- 
dition to this requirement, fixed points are obtained when 
TV = [V]. So to investigate the conditions needed for de- 
terming the fixed points of the system, the second time deriva- 
tive of the function, N = N (t), is obtained and set down in a 
non-autonomous form as 



N 



1 ~ AXt{X l~^ + 2X 2 (Xt-pf-X 2 

fl + t 1 



A. (4) 



A look at Eqs. (O and reveals the existence of a fixed 
point at A = N* = 0, when the conditions N = A = 
are simultaneously satisfied. This fixed point is obtained, as 
Eq. (O indicates, when t — > ±00. The existence of other 
fixed points will depend on whether or not Eqs. (|3) and © 
will have common roots in t when N — N = 0. Solving for 
t under these conditions will also shed much light on how // 
and j3 influence the dynamics of the bimodal system. 

The most direct way to appreciate the role that /1 plays is 
to study it for the case of (3 = 0. In this perfectly symmetric 
limit, on solving the explicitly time-dependent part of N = 0, 

three roots of t are obtained, one at t = and the other two 

1/2 

at t = ± (A~ 2 — fij . It shall be very important to note at 
this point that none of these three roots of t leads to N = 0. 
The only conclusion that can, therefore, be drawn from this 
fact is that there is only one fixed point in the bimodal sys- 
tem, at N = N* = 0. All other solutions of = give 
only turning points in the bimodal system. For t = one ob- 
tains a turning point of N at No = Afj,, as it can be seen from 

1 /2 

Eq. Q2). The two non-zero roots of t at t = ± (A~ 2 — /1) , 
will give two coinciding turning points on the line, N = 0, 
in the phase portrait. From Eq. (|2), these two turning points 
are to be found at At = AA~ 2 exp (/iA~ 2 — l). Now it has 
already been argued in Section|Ill]that to have three real roots 
of t, it is necessary to fulfil the condition, fi < A~ 2 . By means 
of this condition it can also be argued that At > No. Hence, 
as long as /i < A~ 2 , there shall remain three turning points, 
one at N = Nq and the other two coinciding at N = At- As 
fi — > A~ 2 , these two turning point positions approach each 
other in the N—N phase portrait. When fi = A~ 2 , there is a 
bifurcation-like behaviour because of which the turning point 
at Ao and one of the turning points at At annihilate each 
other. What survives this mutual annihilation is the second 
"concealed" turning point at At- Subsequently, for /1 > A~ 2 , 
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FIG. 3: A symmetric phase portrait (with ft — 0), using the param- 
eter values of A and A taken from Fig. [TJ When fi — 0, the phase 
trajectories are shown by the continuous curve at the extreme left. 
The closed loop has the largest area in this case. When [i — 0.5A -2 , 
the loop size shrinks, as the dotted curve shows in the centre of the 
plot. When ji = A -2 , the loop vanishes as the curve on the extreme 
right shows. The existence of a closed loop in the plots above is a 
signature of bimodality. When fi = A~ 2 , a bimodal-to-unimodal 
transition occurs, and the closed loop disappears, i.e. bimodality is 
lost. All along there is a fixed point at N = N = 0, joined to itself 
by a homoclinic trajectory. 



only this turning point will continue its existence, and the dis- 
tribution will become unimodal. 

The entire sequence of this behaviour of the turning points 
over the range, < /i < A~ 2 , has been depicted in Fig. [3] It 
is obvious from this plot that bimodality is intimately related 
to the existence of a closed loop for a phase trajectory in the 
N—N portrait. This is one clear message that has emerged 
from the phase-portrait analysis. When /j, — > 0, bimodality 
becomes very pronounced, and when fi — > A~ 2 , bimodal- 
ity becomes progressively enfeebled, disappearing altogether 
when [i = A~ 2 . It is to be noted that in the former case, the 
sole fixed point at N = N* = coincides with the turning 
point, No = 0, while in the latter case, bimodality is erad- 
icated because of a mutual annihilation between the turning 
points, No 0) and A^t- Once bimodality disappears and 
the distribution becomes unimodal, another fact that stands 
out clearly through the phase-portrait analysis is the existence 
of a homoclinic trajectory that joins the fixed point, N* = 0, 
to itself. A relevant matter to consider here also is that the 
closed loop that traces bimodality in the phase portrait con- 
veys the impression of a local periodic behaviour (which ef- 
fectively implies that the two peaks in the bimodal distribution 
give a local scale), and when /i — > 0, this local periodicity of 
what is now a strongly bimodal system, takes place on a time 
scale of A -1 . 

Equipped with an understanding of how /j, influences the 
dynamics of a bimodal system, the next effect to consider is 
the one of breaking the symmetry that is maintained for (3 — 0. 
Apropos of this, the first thing to note is that when /3 / 0, 
the two hitherto coinciding turning points at N = Nt are 



separated from each other. A simple perturbative argument 
suffices to demonstrate this. If /3 1, then about the turning 
point, A^t (determined for t 2 = A~ 2 — fi), it is possible to 
linearise all terms involving j3 in Eq. (0, and obtain the two 

r 1 /2~ 

separated roots as A" T ± ~ At 1 ± 2/3 (l - //A 2 ) . 

As j3 increases in magnitude, the position of A^t+ shifts to 
higher values of A" on the line N = {), while the position of 
iVx- slides in the opposite direction. This pair of roots cor- 
responds to the two asymmetric peaks in the bimodal curve, 
with A^t+ tracing the locus of the higher peak, and Nt- the 
lower. Between the two of them, these peaks flank a point of 
local minimum in the bimodal distribution, which itself turns 
out to be the third turning point on the line, A" = 0. When 
,3 = 0, this point is at A" = N , with A^ < At, but when 
,3/0, the position of this third unpaired turning point shifts 
to larger values of N, and for sufficiently high values of j3, 
it suffers a separate pair annihilation with the shifting turn- 
ing point at A^t- ■ The argument in favour of this contention 
comes from solving for A" = in Eq. ||3), specifically the 
non-autonomous, time-dependent part. What transpires is a 
third-order polynomial in the variable, t, which goes as 
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(5) 



Performing a variable transformation, t = r+h, in the forego- 
ing relation, and setting h = f3 (3A) , one obtains the canon- 
ical form of a cubic equation, 

r 3 + Pt + Q = , (6) 

whose analytical roots are extracted by the application of the 
Cardano-Tartaglia-del Ferro method. Noting that 
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in the present case, the discriminant of the cubic equation can 
be expressed as 



4 



p3 

27 



(7) 



The sign of T> is crucial here. If T> > 0, there will be 
only one real root of r (the other two being complex, which 
are known to occur in pairs). On the other hand, if T> < 0, 
all the three roots of r will be real. Now, knowing t, leads 
easily to knowing t, which in turn, on using Eq. ©, gives the 
corresponding value of N. This will be a turning point on the 
line N = 0. If there are three real roots of r (for V < 0), then 
there will also be three turning points of A^ (t). If, however, 
there is only one real root of r (for T> > 0), then there will 
be only one turning point of A^ (t), the other two having been 
annihilated on the line, A^ = 0. 

Evidently, a critical condition is reached when V = 0, at 
which a bimodal-to-unimodal transition occurs. This might be 
viewed as a bifurcation in the bimodal system, for the value 
of j3 — /3b. To understand how asymmetry in the bimodal dis- 
tribution brings about this type of bifurcation behaviour, it is 
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FIG. 4: When symmetry is broken (for f} ^ 0), there are three 
distinct turning points for N = 0. When the asymmetry is large, two 
of these fixed points are lost due to a bimodal-to-unimodal transition. 
This happens when the discriminant of Eq. ((6), T> = 0. The variation 
of T> against |/3| (since /3 < for the westward traffic) shows that 
the transition occurs when |/3| = |/?b| — 0.21 for the eastward traffic 
(traced by the dotted curve), and when \/3\ = \/3b\ — 0.46 for the 
westward traffic (the continuous curve). For both of these cases the 
values of n and A are what they are for Figs. [T]&|2] respectively. 



a simple task to follow how V depends on /3. This variation 
has been plotted in Fig.Hl using the values of // and A needed 
to calibrate the data of both the east and west-moving traffic. 
Considering |/3| (since j3 < for the west-going traffic), a 
transition from a bimodal distribution to a unimodal distribu- 
tion (or vice versa) can be seen to take place (when T> = 0) 
for | /3b | — 0.21 in the case of the eastward traffic, and for 
|/3b| ~ 0.46 in the case of the westward traffic. 

Interestingly, going back to Figs. Q~|&[2] one can see that 
in the former case, the value of |/3| needed to fit the data 
(|/3| = 0.09) is much less than what can cause the bimodal- 
unimodal bifurcation (|/3b| ~ 0.46), while in the latter case 
the reverse is true (here |/3 b | ~ 0.21 < |/3| = 0.24). The con- 
sequence of this is that, as opposed to Fig.Q] there is no con- 
spicuous secondary peak in the bimodal distribution shown in 
Fig. |2] The implications of both these conditions, as far as the 
dynamics on the N—N phase portrait is concerned, have been 
shown in Figs.[5]&|6] These two plots show the dynamics of 
bimodality for the west-moving and the east-moving traffic, 
respectively. In both the cases, when /3 = 0, there is a point on 
the line N — 0, where the curve loops over itself. The shifting 
of this point away from the aforesaid line, suggests breaking 
of symmetry and asymmetric bimodality. For /3 < (west- 
ward traffic), this point shifts below the line, while for /3 > 
(eastward traffic), the point shifts above. So to give a suc- 
cinct description, it can be said that fi causes bimodality, and 
/3 causes a breaking of symmetry. Taking these two aspects 
of the two parameters together, what results is an asymmetric 
bimodality. And changing the values of both these parameters 
can cause a bimodal-to-unimodal transition in the distribution, 
much in the likeness of a bifurcation. 
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FIG. 5: Phase plot of the bimodal distribution of the traffic flow due 
west. Both the curves are plotted for the values of A, n and A used 
in Fig.Q] When /? = 0, the resulting plot is shown by the continuous 
symmetric curve. When j3 = —0.09 (as in Fig. [I), asymmetry sets 
in, as shown by the dotted curve. 
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FIG. 6: Phase plot of the bimodal distribution of the traffic flow 
due east. The plotting here follows the scheme of Fig. [5] Both the 
curves are plotted for the values of A, /i and A used in Fig. [2] The 
continuous curve shows the symmetric case (for /3 = 0), while the 
asymmetric dotted curve is plotted for = 0.24 (as in Fig.©- bince 
this value of ft > |/3b | (the latter value is to be read from Fig. |4}, a 
bimodal-to-unimodal transition has already taken place, and instead 
of a closed loop on the dotted curve, all that remains is a cusp just 
above the line, N = 0. 



V. CONCLUDING REMARKS 

An explicit model equation for bimodal behaviour notwith- 
standing, the dynamics and the phase-portrait analysis have 
helped further in unearthing some general properties of a bi- 
modal distribution. The first is that the signature of bimodal- 
ity in a phase-portrait is the existence of a closed loop, and the 
second is the existence of a homoclinic trajectory connecting 
a lone fixed point to itself in the phase portrait. One might 
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suppose that these are features that will hold true in any other 
kind of bimodal distribution. 

One can make the case that modelling of vehicular traffic 
flows in the manner done here can help in designing smarter 
traffic control mechanisms, and afford a new perspective on 
addressing traffic-related problems — as, for instance, adjust- 
ing and allocating an appropriate amount of green-signal time 
when traffic is congested along one lane, and moving freely 
along another. 

The model presented here is based on a relatively uncom- 
plicated bi-directional flow of traffic along one axis only (east 
to west and vice versa). There is scope for refinement, when 
one considers the complications associated with the vast traf- 
fic networks of very big cities. With the aid of advanced mod- 



elling techniques, it should become possible to simulate and 
analyse conditions for which data may not be available, or for 
situations in which it may be difficult to perform realistic ex- 
periments on large scales. All of these would go a long way 
in building a generalised model that could be applied to traffic 
flows in any arbitrary direction. 
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